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ABSTRACT 


A  preliminary  coldroom  experiment  was  conducted  to  determine  the  diffusion  rate  of  several  reactive 
atmospheric  chemical  contaminants  in  snow  at  -7°  C.  These  species  included  formaldehyde,  p-nitrotoluene 
(PNT),  HNO3,  and  HCl.  Because  of  detection  or  contamination  problems  we  could  not  quantify  diffusion 
rate  for  formaldehyde,  PNT,  or  HNO3.  However,  the  diffusivity  of  HCl  in  snow  (p  =  0.38  ^cm^)  was 
found  to  be  >  7.7  X  10^  cm^  /s.  This  is  about  two  orders  of  magnitude  greater  than  hitherto  expected  and 
suggests  tiiat  HCl  may  be  rapidly  redistributed  or  lost  from  seasonal  snowpacks.  Thus,  one  should  be 
cautious  in  interpreting  temporal  records  of  HCl  and  other  volatile  species  in  seasonal  or  polar  snow- 
packs,  until  their  integrity  is  better  established. 
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A  Preliminary  Investigation  of  Diffusion 
of  Volatile  Atmospheric  Acids  in  Snow 
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1  INTRODUCTION 

Trace  constituents  of  the  atmosphere  can  contribute  to  important  processes 
such  as  ozone  destruction,  acid  precipitation,  and  global  warming.  It  is  important 
to  understand  the  sources  and  fates  of  these  compoimds  in  order  to  assess  the  role 
of  man’s  activities  on  these  processes.  Some  trace  species  that  have  been  impli¬ 
cated  are  the  mineral  acids  hydrochloric,  sulfuric,  nitric,  and  low  molecular 
weight  oxygenated  species,  formate,  acetate,  formaldehyde,  sulfur  dioxide,  and 
hydrogen  peroxide.  These  species  are  very  soluble  in  water  and  are  efficiently 
scavenged  and  dissolved  by  small  water  droplets  in  the  atmosphere.  When  these 
droplets  are  quickly  frozen,  as  often  occurs  in  the  atmosphere,  the  trace  species 
are  “frozen  m”  (Iribame  and  Pyshnov  1990)  at  concentrations  exceeding  their 
equilibrium  solubility  in  ice. 

Permanently  frozen  precipitation  may  then  serve  as  a  semi-permanent  record 
of  historic  atmospheres,  and  polar  ice  core  analysis  could  be  used  retrospectively 
to  establish  past  conditions.  One  may  doubt  fliis  conjecture,  however,  as  volatile 
species  would  be  expected  to  redistribute  over  time  so  as  to  achieve  local  ice/air 
equilibria.  We  would  expect  these  changes  to  be  reflected  in  temporal  changes 
in  volatile  species  concentrations  in  snowpacks,  and,  in  fact,  these  kinds  of 
changes  have  been  observed  in  mid-latitude  snowpacks  (Hogan  et  al.  1985) 
and  in  laboratory-stored  snow  (Cragin  and  McGilvary  1995).  The  kinetics  of 
redistribution  processes  and  the  ultimate  fate  of  various  contaminants  remain 
obscure,  however. 

Contaminants  also  enter  snowpacks  by  diy  deposition  on  the  fallen  snow. 
This  applies  not  only  to  atmospheric  sources  but  to  sources  such  as  soil  beneath 
the  snow.  Because  of  the  prevailing  temperature  gradients,  volatile  contaminants 
tend  to  accumulate  in  the  coldest  snow,  which  is  often  the  top  of  the  pack  (Hogan 
et  al.  1985,  Hogan  and  Leggett  1995).  Because  these  phenomena  are  hard  to 
study  in  nature,  most  of  our  information  comes  fi'om  laboratoiy  investigations. 
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However,  no  direct  information  on  diffusion  of  trace  volatile  species  in  snow 
was  found.  Diffusion  of  solute  vapors  through  a  porous  medium  such  as  snow  is 
expected  to  be  retarded  with  respect  to  vapor  movement  in  air  by  interaction  with 
the  media  and  tortuosity  effects. 

Available  information  suggests  that  interaction  of  acidic  species  with  ice  par¬ 
ticles  is  strong,  and  that  diffiision  into  the  bulk  ice  occurs.  The  diffusivities  of 
hydrochloric  and  nitric  acids  from  dilute  v^or  sources  in  ice  have  recently  been 
estimated.  The  reported  values  range  from  about  10"'^  cm  Vs  at  -8°  C  for  ice 
single  crystals  (Domine  and  Thibert  1995)  to  aroimd  cmVs  for  polycrys¬ 
talline  ice  (Diehl  et  al.  1995).  The  latter  value  would  presumably  be  more 
relevant  for  atmospheric  studies  and  for  the  work  reported  here  since  snow 
is  essentially  unconsolidated  polycrystalline  ice.  If  so,  we  would  anticipate 
essentially  complete  permeation  of  the  ice  spheres  during  the  time  course  of 
our  experiment,  although  complete  equilibration  with  the  ice  lattice  could  take 
years,  assuming  the  lower  diffusivity. 

To  avoid  complicating  effects  of  temperature  gradients,  we  used  a  single 
isothermal  temperature  for  this  preliminary  study  of  dilute  vapor  migration 
through  snow.  Compormds  initially  selected  for  study  were  HCl,  HNO3,  for¬ 
maldehyde,  and  p-nitrotoluene.  The  latter  was  added  as  a  reference  compound 
as  one  of  us  had  successfully  used  it  before  in  outdoor,  temperature  gradient 
experiments  (Hogan  and  Leggett  1995). 
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2  EXPERIMENTAL  METHODS 

All  laboratory  work  was  performed  in  a  Class  100  clearu’oom  or  Class  100 
clean  air  station  in  a  coldroom.  Because  of  the  finite  background  concentrations 
of  many  atmospheric  species  in  natural  snow,  we  prepared  clean  artificial  snow 
from  frozen  Milli-Q  deionized  water  as  described  below. 

First,  we  added  Milli-Q  water  to  pre-cleaned  7.5-cm-diameter  polyethylene 
tubes  that  were  heat-sealed  on  the  bottom.  Leaving  about  20  cm  of  air  space  for 
expansion  above  the  60  cm  of  water,  we  sealed  the  top  and  slowly  (1.8  mm/hr) 
froze  the  water  vertically  by  lowering  it  into  a  circular  opening  in  a  freezer  (see 
Cragin  and  McGilvary  [1995]  for  a  more  detailed  description  of  the  apparatus 
and  process).  We  sliced  the  resulting  ice  core  into  1-cm-thick  sections  and  rinsed 
them  in  Milli-Q  water  to  remove  surface  contaminants  and  then  mechanically 
crushed  the  sections  in  a  pre-cleaned  Plexiglas  chamber  (similar  to  a  mortar  and 
pestle). 

We  sieved  the  resulting  snow  so  that  the  0.425-  to  0.850-mm-size  range 
fell  natmally  through  a  funnel  into  a  10-cm-diameter  Plexiglas  column.  We 
suspended  separate  1  -mL  vials  of  formaldehyde,  0. 1  M  HCl,  0. 1  M  HNOs,  and 
PNT  in  the  airspace  above  the  snow  and  finally  sealed  the  column  with  an 
airtight  cap.  Vapor  was  free  to  move  within  the  chamber  from  the  xmcapped 
source  vials  to  the  snow.  We  left  the  apparatus  in  a  coldroom  at  -7°  C. 

After  one  week  we  opened  the  column,  removed  the  source  vials,  and 
sectioned  the  snow  core  crosswise  in  a  clean  air  station  in  the  coldroom.  The 
resulting  1-  or  2-cm-thick  wafers  were  melted  and  analyzed.  Analyses  for  for¬ 
maldehyde  and  PNT  were  done  using  gas  chromatography  while  chloride  and 
nitrate  were  determined  using  ion  chromatography. 
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3  RESULTS  AND  DISCUSSION 

We  failed  to  detect  p-nitrotoluene  (PNT)  in  any  of  the  snow  samples.  This 
was  somewhat  surprising  because  it  had  been  detected  in  earlier  work  (Hogan 
and  Leggett  1995)  in  snow  over  PNT-spiked  soil  at  a  similar  temperature  and 
using  the  same  method  of  analysis.  The  difference  was  that  the  earlier  experiment 
was  conducted  outside  in  the  presence  of  a  strong  temperature  gradient,  so  that 
the  source  was  at  a  higher  temperature  (~0°  C)  compared  to  -7°  C  in  the  present 
experiment.  The  lower  temperature  is  expected  to  reduce  the  somce  strength 
by  a  factor  of  between  2  and  3.  The  temperature  of  the  snow  in  which  PNT  was 
detected  in  the  former  experiment  may  also  have  been  a  bit  colder  than  -7°  C, 
and  we  failed  to  detect  it  when  the  snow  was  warmer  dian  -5°  C.  Evidently,  the 
sensitive  temperature  dependence  of  PNT  source  strength  and  retention  slow  its 
diffusion  under  isothermal  conditions.  This  seems  likely  to  be  a  very  general 
phenomenon  rather  than  unique  to  PNT. 

We  also  did  not  detect  nitrate  in  any  of  the  snow  samples,  but  this  may  be 
because  concentrations  were  below  the  instrumental  detection  limit.  Because  of 
its  high  vapor  pressure,  we  expected  nitric  acid  to  volatilize  and  diffuse  readily 
through  the  snow,  but  if  nitrate  concentrations  were  below  20  pg/L,  they  would 
have  been  undetected. 

Attempts  to  determine  formaldehyde  likewise  proved  frustrating.  The  prob¬ 
lem  was  not  in  failure  to  detect,  however,  but  in  imexpected  contamination. 
Unfortunately,  the  original  ice  used  to  make  the  snow  for  the  column,  which  was 
not  analyzed  until  after  the  experiment,  was  contaminated.  Similar  amounts  were 
also  found  in  all  column  samples,  which  precludes  drawing  any  conclusion  ex¬ 
cept  that  the  results  are  not  inconsistent  with  formaldehyde  diffusion  through 
the  column.  A  definitive  result  must  await  further  experimentation.  Analysis  of 
several  other  batches  of  lab  Milli-Q-purified  water  also  had  varying  amoimts  of 
contamination.  One  sample  of  well  water  fi-om  a  local  source  was  analyzed  and 
revealed  undetectable  levels  of  formaldehyde. 

Before  future  experiments  are  conducted,  potential  sources  of  laboratory 
contamination  other  than  water,  such  as  exposure  to  contaminated  lab  air,  must 
also  be  ruled  out.  Indoor  air  generally  contains  more  formaldehyde  than  outdoors 
due  to  degassing  of  materials  containing  formaldehyde  in  a  confined  space.  Since 
formaldehyde  occurs  ubiquitously  in  the  environment  it  will  probably  prove  im¬ 
possible  to  completely  eliminate  diis  source  of  contamination;  however,  it  should 
be  possible  to  do  the  experiment  by  minimizing  e?qposure  to  sources  of  contami¬ 
nation  and,  if  necessary,  by  increasing  the  formaldehyde  source  strength. 
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Figure  1  shows  the  concentration  profile  for  HCl  at  the  end  of  seven  days. 
The  profile  is  atypical  in  that,  although  the  gas  has  permeated  through  the  entire 
column,  the  concentrations  behind  the  fi'ont  decrease  regularly.  With  a  constant 
source,  uniform  concentrations  behind  the  front  would  be  expected.  The  simplest 
explanation  for  this  result  is  that  the  source  did  not  remain  constant  over  time  but 
became  depleted.  This  is  prone  to  h^pen  with  unstirred  liquid  solution  sources 
such  as  used  here,  and  would  be  exacerbated  by  the  restricted  geometry  of  the 
source,  because  of  the  mismatch  between  the  rates  of  gas  phase  vs.  liquid  phase 
diffusion.  Diffusion  coefficients  of  solutes  are  roughly  four  orders  of  magnitude 
higher  in  air  than  in  aqueous  solutions.  This  can  cause  depletion  of  the  surface 
l^er  of  the  source  and  proportionately  reduced  solute  vapor  flux  when  diffusion 
through  the  porous  medium  is  as  rapid  as  it  was  here.  Rapid  diffusion  of  HCl 
through  the  snow  coliunn  was  not  anticipated,  as  the  literature  suggested  re¬ 
tention  by  the  ice  surface  would  be  more  efficient  (Diehl  et  al.  1995,  Domine 
and  Thibert  1995). 


HCl  Diffusion  in  Snow 


Figure  1.  Chloride  concentration  in  snow  versus  distance  from  the  source 
(HCl).  The  total  length  of  the  snow  column  was  30  cm. 

One  might  ask  why  the  HCl  doesn’t  diffuse  back  through  the  snow  column. 

It  does.  Because  diffusion  occurs  in  both  directions  (it  actually  occurs  in  three 
dimensions,  but  the  column  wall  restricts  radial  movement),  HCl  must  have 
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diffused  backward  out  of  the  snow  once  the  source  became  depleted.  The  lower 
concentrations  closer  to  the  source  probably  reflect  this  as  the  HCl  diffused 
out  of  the  snow  and  possibly  repartitioned  with  solutions  in  vials  of  the  other 
sources.  Once  the  HCl  became  depleted  in  its  somce  vial,  the  solutions  in  the 
other  vials  became  sinks  for  HCl  in  the  headspace  and  the  snow. 

Since  the  diffiision  front  passed  completely  through  the  column  we  cannot 
calculate  an  exact  diffiisivity  of  HCl  for  the  experiment.  We  can  estimate  a  mini¬ 
mum  value,  however,  assuming  the  following  equation  applies; 


r  =  (2£>0‘^^  (1) 

where  r  is  the  distance  of  the  front  from  the  source  at  the  end  of  the  e3q)eriment, 
D  is  the  effective  diffiisivity  or  diffusion  coefficient,  and  t  is  time.  In  this  case 
r  >  30  cm  and  t  =  5.83  x  10^  s  (6.75  days).  Solving  for  D  we  obtain  a  value  of 
/)  >  7.7  X  10"^  cm^/s.  This  calculation  shows  why  the  HCl  flux  through  the  col¬ 
umn  was  limited  by  mass  transfer  in  the  source  solution.  Diffiisivity  of  HCl  in 
water  at  -7°  C  is  expected  to  be  no  more  than  10"^  cmVs  (Thibodeaux  1979). 
Since  this  is  about  two  orders  of  magnitude  less  than  its  minimum  estimated 
diffiisivity  in  the  snow  column,  it  is  easy  to  see  why  such  a  mass  transfer 
limitation  occurred. 

Diffusion  of  gases  in  porous  media  is  retarded  relative  to  their  diffusion  in 
air  by  sorption  on  the  media  and  tortuosity  effects.  We  estimate  the  effect  of 
tortuosity  by  assuming  spherical  geometry  of  the  snow  particles  and  additivity 
with  respect  to  the  two  volume  fraetions,  the  air-filled  voids,  and  the  particle 
fraction.  If  we  denote  these  by  Vy  and  Ps,  the  tortuosity  factor  would  be  just 
Vy  +  Ps  (pi/2).  The  measured  bulk  density  of  the  snow  was  0.38  g/cm^.  Therefore 
an  estimate  of  tortuosity  for  the  experiment  is  0.62  +  0.3  8  (pi/2)  =  1 .22.  This  is 
probably  a  minimum  estimate  due  to  the  crude  assumptions,  but  seems  in  reason¬ 
able  agreement  widr  values  reported  for  other  porous  media  (Thibodeaux  1979). 

It  was  of  interest  to  estimate  a  sorption  coefficient  for  HCl  on  snow.  In  prin¬ 
ciple  the  bulk  sorption  coefficient  can  be  determined  from  the  effective  diffiisivi¬ 
ty  and  tortuosity.  We  can’t  do  this  since  we  don’t  have  a  definite  value  for  the 
diffiisivity,  but  we  can  estimate  a  maximum  value  by  assuming  that  r  =  30  cm, 
i.e.,  the  length  of  the  snow  column.  The  molecular  diffiisivity  of  HCl  in  air  at 
-7°  C  was  estimated  as  0.14  cm  Vs  from  reported  values  for  other  compounds 
of  similar  molecular  weight  at  0°  C  (Thibodeaux  1979;  the  lower  temperatme 
should  have  negligible  effect).  We  divide  this  by  the  tortuosity  to  obtain  a 
theoretical  diffiisivity  for  HCl  in  the  snow  column  in  the  absence  of  adsorption. 
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Our  estimated  unretarded  diffusivity  is  then  about  0.1 1  cmVs.  An  equation  for 
retarded  diffusion  in  porous  media  is 

R^DJD,  =  \+{pj:iJV,)  (2) 

where  is  the  effective  diffiisivity  without  retardation,  A  is  the  retarded  dif¬ 
fusivity,  >  7.7  X  cm  Vs,  Ps  is  now  the  bulk  density  of  the  snow  in  g/cm^,  Fy 
is  the  void  fraction  in  mL/cm^,  and  Kb  is  the  bulk  sorption  coefficient  and  has 
dimensions  of  mL/g.  Rearranging  in  terms  of  Ab  we  have 

Kb  =  {R-\)V,lp,.  (3) 

Substituting  the  above  values  gives  Kb  <  140  mL/g. 

An  alternative  expression  for  sorption  coefficient  is  in  terms  of  surface  area. 
To  do  this  we  just  divide  Kb  by  the  surface  area  of  the  snow  particles  in  cmVg. 

We  don’t  have  precise  knowledge  of  this  but  can  estimate  it  from  the  sieve 
sizes  used  to  prepare  the  snow  and  assiuning  spherical  geometry  of  the  particles. 
For  the  sieve  openings  we  used  the  particles  range  from  0.425  to  0.850  mm  in 
diameter.  This  gives  a  range  of  specific  siuface  areas  of  79  to  158  cmVg.  For 
estimation  purposes  we  use  the  mean  or  120  cmVg  and  estimate  a  surface  area- 
normalized  sorption  coefficient  as  A  <  1-2  cm  (actually  cm^  air/cm^  ice).  This 
is  the  more  meaningful  measure  if  HCl  is  adsorbed  on  the  ice  surface  rather  than 
dissolved  in  the  bulk. 

Actually  the  maximum  concentration  measiu'ed  in  the  current  experiment 
of  approximately  16  ng/mL  represents  only  1.3  x  10"®  g/m^,  which  is  at  most 
1%  of  monolayer  coverage.  Thus  we  have  no  reason  to  infer  sorption  into  the 
bulk  or  postulate  uptake  into  a  liquid-like  surface  l^er,  in  agreement  with  others 
(Domine  et  al.  1995). 

Of  interest,  and  as  a  check  on  the  reasonableness  of  the  above  calculation, 
we  can  back-calculate  the  minimum  HCl  vapor  concentration  that  should  have 
been  obtained  in  the  column.  Since  Kb  =  CyJCa  =140  mL/g,  where  Cb  and  Ca  are 
the  respective  HCl  concentrations  in  bulk  ice  and  air,  we  can  solve  for  Ca.  We  get 
Ca=  1.1  X  10"®  g/mL  or,  in  terms  of  pressiu’e,  5.0  x  10"®  torr.  This  appears  to  be 
too  high  based  on  available  experimental  data  and  recent  fitting  models  (Miller 
1983,  Luo  et  al.  1995).  The  model  supports  a  partial  pressure  for  HCl  in  equili¬ 
brium  with  3.7  wt%  HCl  in  water  at  -7°  C  of  ~1.0  x  10"®  torr. 

However,  if  the  model  value  is  correct,  Kb  becomes  7.3  x  10®  mL/g  and  D, 

2.5  X  10"®  cmVs.  At  this  diffusivity  HCl  should  have  permeated  only  to  about  5 
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cm  in  the  column  instead  of  all  the  way  through.  Resolution  of  this  discrepancy 
must  await  future  experiments.  It  is  conceivable  that  the  other  species,  namely 
HNO3  and  formaldehyde,  interfered  with  sorption  of  HCl  on  the  snow.  We  note 
that  the  partial  pressure  of  HNO3  under  the  experimental  conditions  (estimated 
according  to  Luo  et  al.  1995)  would  have  been  nearly  an  order  of  magnitude 
higher  than  that  of  HCl,  so  it  may  have  competed  for  sorption  sites  on  the  ice 
surface.  We  anticipate  follow-up  experiments  in  which  only  one  species  will  be 
present  at  a  time,  and  will  attempt  to  measure  its  actual  partial  pressure  experi¬ 
mentally.  This  is  especially  important  for  these  strong  acids  in  water  as  their 
activities,  i.e.,  partial  pressures,  are  highly  concentration-dependent. 
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4  ENVIRONMENTAL  IMPLICATIONS 

Although  this  work  is  preliminary,  the  rapid  diffusion  of  HCl  observed  in 
snow  at  -7°  C  suggests  that  this  and  possibly  oflxer  mineral  acids  may  be  rapidly 
redistributed  and  lost  from  seasonal  snowpacks.  This  might  also  explain  why 
some  acid  anions  are  lost  on  storage  of  snow  samples  (Cr^in  and  McGtilvary 
1995).  In  nature,  temperature  gradients  tend  to  drive  contaminants  upward 
through  snowpacks  (Hogan  and  Leggett  1995)  where  concentration  gradients 
can  return  them  to  the  atmosphere  by  diffusion  and  convection.  Thus,  one 
should  be  cautious  in  interpreting  measurements  of  potentially  volatile  species 
in  seasonal  snows  or  in  polar  fim,  until  their  integrity  is  better  established. 
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